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Abstract 

The up-dated three-phase concept [jy] of nuclear matter evolution in course of 
cooling down - from the phase of quark-qluon plasma (QGP) through the interme- 
diate phase allowing for massive constituent quarks Q (valons), pions and kaons 
(Q7rK) to the phase of hadronic matter (H) - is exploited for the treatment of rela- 
tive hadronic yields in the central region of heavy ion collisions. The most attention 
is paid to the description of the Q-7rK-phase which is argued to be a gaseous one and 
lasts until the valonic spacing approaches the confinement radius (at the tempera- 
ture Th — (110 ± 5) MeV), when the valons start fusing to be locked, in the end, 
within the hadrons. The hadronic yields emerged from thermal treatment of Q-7rK- 
phase and simple combinatorial approach to the hadronization process are shown 
to fit the available experimental data from AGS, SPS and RHIC quite well. This 
provides an alternative insight into the real origin of the observed relative hadronic 
yields which is (to a considerable extent) free of the well known puzzle inherent in 
some conventional models where the early chemical freeze-out is assumed: namely, 
why the gaseous thermal approach to actually tightly packed (even overlapping) 
hadrons seems workable? Many predictions for the other hadronic yields which 
could be observed at these machines as well as at LHC are given. 
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Introduction 

The conventional way of treatment of hadron production in the central rapidity region 
of heavy ion collisions at high energies incorporates a QGP as the short initial phase 
of nuclear matter (just after the nuclei encounter each other) which then undergoes the 
chiral breaking transition into hot and dense hadronic matter, the latter evolving somehow 
into experimentally observed free hadrons. The relative yields of different species of 
these hadrons depend on the equation of state (EoS) of the nuclear matter they are 
originated from. It has been noticed || that the SPS data on the hadronic yields are 
quite compatible with the assumption about the early chemical freeze-out in an ideal gas 
of crucially modified in-medium hadrons (their effective radii has being assumed to be as 
small as ~ 0, 3 fm only!) at the temperature T c h ~ 170 MeV [] (which is rather close to 
the chiral breaking temperature itself, T c ~ 200 MeV). The theoretical observation has 
been also made B that thermal description of such a type could be spread to elucidate the 
AGS and SIS data as well under the assumption that chemical freeze-out temperature T c ^ 
was varied properly along with baryonic chemical potential of fireball (expanding nuclear 
medium), being about 110 MeV and 60 MeV for AGS and SIS, respectively. As a result, 
a suggestion about possible existence of low temperature QGP has been put forward and 
then discussed in some more detail || ||. Both attributes of this approach - crucially 
modified in-medium hadrons (but still survived!) and low temperature chiral transition - 
seem rather controversial and mysterious. That is why a feeling of dissatisfaction remains 
and looking for some more consistent approach seems quite reasonable. 

In this connection, the attention seems to be called reasonably to the rather old- 
fashioned notion of valons Q linked properly to the corresponding current quarks - first 
of all, to the valons Q(q), q = u,d (mQt U) d) = m Q{ q ) — 330 MeV), and Q(s) {mQ( s ) — 
480 MeV) - which was widely and fruitfully exploited before QCD was developed as a 
consistent field theory that had no valons as its inalienable entities. However, till now, 
QCD remains actually a quite workable theory of hard processes, but it suffers from the 
lack of even qualitative results for the processes at low and intermediate energies, i.e., just 
for those ones which proceed within nuclear matter below the chiral breaking temperature 
T c . Thus, maybe, embedding the valons unambiguously into the body of QCD is still 

^he straightforward estimate shows that no nearly ideal gas of slightly modified hadrons could ex- 
ist at the relevant particle (hadron) densities, since the normal hadron wave functions would overlap 
substantially. 
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waiting for its turn, although the known endeavors met no success |J. Qualitative physical 
motivation in favor of possibly essential role of valonic mass scale and of an intermediate 
phase, allowing for valons and separating QGP and hadronic ones, can be found in ref.s 
|6|, [7], |8[]. It is worthy to mention that experimentally observed (in AA collisions) significant 
excess in low-mass dilepton (e + e~) yield (as compared to pA ones) could be described 
quite well || in the framework of an approach allowing for Q7rK-phase. 

Recently the general idea has been put forward |IJ which could provide an insight 
into the physical reasons underlying the unified thermal description of yields of different 
hadron species observed in heavy ion collisions. In the present paper, we show in more 
detail how could valons help in treatment of relative hadronic rates, getting rid of the 
above curiosity inherent in the conventional approaches. The theory is confronted to the 
available results of AGS, SPS and RHIC, and many other ratios of hadronic yields are 
predicted (for LHC as well). 

General Description of the Approach 

Most probably, just after chiral symmetry breaking, the nuclear matter is subjected 
for dual description ||: it can be treated either as highly compressed "hadronic liquid" or 
as a state, in which the most of hadron species can not survive and the dominant degrees 
of freedom are not still hadrons but valons instead Q. Of course, one can choose either. 
However, for taking the first way, one is supposed to know, at least, the corresponding EoS. 
Who is aware of it? The basic point of the suggested approach is that the valons which 
would be produced in course of chiral symmetry breaking are to be treated reasonably as 
a gas from the very beginning (in contrast to the hadrons!) due to their it really small size 



nfl, r ~ 0,3 fm. Indeed, even at T = 170 MeV the particle density within the "ideal" 



Boltzman valonic gas would be about 

^ [2( ^ fK2( ^m )msh{llmlT) + ( ^ fK2 (^ )coshMT)] s 16n -3 

where fiQ( q ) and fiQ( s ) are typical for heavy ion collisions Q(q)- and Q(s)-valon chemical 
potentials, respectively^, and thus the "valonic bodies" themselves occupy about 10% of 
the total volume only. This rather optimistic estimate noticeably worsens after taking into 
2 Presence of massive valons is just indicative for chiral symmetry breaking. 

3 We do not distinguish here between chemical potentials of Q(u)- and Q(c?)-valons caused by difference 
in proton and neutron content of heavy colliding nuclei, although at low energies and large values of [i q /T 
(SIS, AGS) it may result in a noticeable effect. 
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account some balancing fraction of the "large-sized" pions and kaons which are produced 
inevitably as the chemical equilibrium sets in [], see Fig. 1. Fortunately, this hadronic 
fraction never exceeds 25% (see below), and thus this three-component Q7rK state is still 
to be treated reasonably as a gas Q, although quasi-ideality of this gas becomes somewhat 
more questionable. That is why below we proceed as far as possible with a more general 
consideration based on the chemical kinetics only, and then confront the results to what 
one can obtain in the selfconsistent ideal gas approximation. We will see a 'posteriory that 
the results emerged from these two ways of doing are quite compatible. 

Now, the general pattern of evolution of the hot nuclear fireball produced in course 
of heavy ion collisions looks as follows. While expanding and cooling down from T c to 
Tjj, the nuclear matter is gradually enriched with pions and kaons and impoverished 
with valons. Situation changes dramatically when mean spacing of valons scales up to 
become about the confinement radius. As a result, gradually increasing color screening 
length approaches the critical value and the dominant degrees of freedom become no 
longer valons (since they are getting not free) but the hadrons themselves, i.e., the bulk of 
hadrons is to build up. This is a phase transition, if the hadronization proceeds near the 
same temperature T ~ Th] otherwise, this is a phase crossover. The lattice calculations 
suggest that peculiarities of phase transition and even its occurrence itself is regulated by 
the value of chemical potential (or net baryonic density) within nuclear matter. We do 
not touch here this delicate question and concentrate on the pattern linked to a certain 
phase transition. This pattern is shown below to face no obvious contradictions with 
the available observables within the nowday accuracy of the data (see the Table). The 
relevant diagrams of hadron production are depicted in Fig. 2. 

Calculations and Results 

At the hadronization phase transition (T = Th), all hadron species are produced, 
generally, in the same manner. Unlike the Q7rK-phase, the correlation length remains 
large comparing to the mean particle spacing (even "infinite", if it is the second order 
one) under the whole this phase transition long. Therefore, multiparticle interactions 

4 A11 other hadronic species possess too large masses and lose their identity within this color conducting 
medium, being very unstable with respect to decay into valons, especially due to collisions; see, e.g., ref. 
for a more consistent treatment of p-meson degradation along with increase of nuclear matter density. 

5 This is especially true near the hadronization temperature Th, what is most relevant in the context 
of the problem under discussion. 
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dominate at this stage. One can assume reasonably also that now, once being created, all 
the hadrons survive because the confinement starts working. 

Thus, the final (observable) pions (kaons) are produced in the two-fold way: first, 
there are equilibrium pions (kaons) mentioned above (those, which are "stored up" in 
course of Q7rK-phase evolution, before overall hadronization, see Fig. 1), and second, 
there are pions coming from Q(q)Q(q) — > 7r + X reactions (kaons coming similarly from 
Q(q)Q{s) and Q{s)Q(q) reactions) just at T = T H . The latter are produced at the stage 
of hadronization only, just like the other hadron species, see Fig. 2. The number of the 
former pions reads obviously 

bn 

where n is the total number of color particles within fireball and b is the pionic fraction 
which is calculated in the Appendix (see eq.s A(l)-A(6)). The numerical value of b 
turns out to be rather stable for SPS, RHIC and LHC: b ~ 0,22, whereas for AGS it is 
considerably lower: b ~ 0,13. Making use of the diagrams shown in Fig. 1(b), 1(c), the 
number of former kaons is estimated quite similarly (see eq. A(7)) to be 

An w + 3n Q(9 ) 

"*+ * * 6(4n, + 3n m ) "«<" (2) 

where uqia is the number of (Q(s) valons. The number of K~ (K°) is obtained by 
replacement here valons antivalons. 

As for the latter pions, one can express their number by tracing a valon Q(q) (antivalon 
Q{q)) within hadronizing nuclear medium. While moving along the mean free path, it 
coalesces with an antivalon (valon) with the probability ^ng^/n (nQ^/n), the total 
rate of such collisions thus being UQ(^nqt q \/n. In each collision, ir- meson plus X are 
produced, "A" allowing possibly for a number j of pions too, see Fig. 2 (the latter is 
also true for the processes of production of other hadrons). Thus, accounting eq. (1), the 
total negative pion yield (which is assumed to be equal to 1/3 of the total rate of pions) 
reads 

N„- K i [-*4 + " + < J »'"W"^ ] + <|Vb + N B + (N K - n K ) + ...] (3) 
o x o n o 



6 Of course, at the hadronization stage, any valonic interaction which could result in producing a 
hadron is reasonably assumed to do produce it. Besides, the valon-hadronic interaction is insignificant 
at this stage. 
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where Nb and Nk are the baryon (antibaryon) and total AT-meson yields, respec- 
tively, and (j) is the mean value of j which is easily estimated to be < j < 1 because of 
the phase space limitation [] (points stand for a number of small terms which are mean- 
ingless within the accuracy of both the theory and the data). Proceeding similarly and 
making use of eq. (2), one obtains for A'-meson yield: 



An* + 3n Q(q) n Q(g) n Q(s) 
(4^ + 3n Q(q) ) Uq{s) In : 

N K o ~ N K o ~ 0,5(^0 + N R0 ), (5) 



Nk + * N K0 ~ - , o x riQ {s) + ^ ™ , (4) 
6(471* + 3nQ (q) ) vv ; 2n 



Nk- and N^o being obtained as before (see eq. (2)). The yields of some other mesons 
are easily estimated by means of the similar combinatorial and/or relevant cross-section 
consideration. For example, 

N ^ n Q{s) n Q (s) 
n 

As for T^-meson, the relative rate N v /N n o can be immediately calculated by comparing 
the corresponding cross sections, since two channels involved {Q{q)Q{q) — > vr7r° and 
Q(q)Q(q) — > vrr/, see Fig. 2(b)) directly compete with each other. In addition to dis- 
tinction in masses of these mesons which results in about (2-=-3)-time contraction (at the 
relevant energies/temperatures) of the 7rr] final phase space (as compared to the nn one), 
one should take into account that no p-meson intermediate state is allowed in 7r ± ?7 chan- 
nel. Being integrated over the energy, the corresponding excess in the tttc yield results in 
an extra factor about 4-f-5, thus making finally N v /N n o ~ (0,08-^-0,1). Unfortunately, the 
flexibility of this estimate is rather large. 

Some additional peculiarities are to be accounted for getting the formulae for the 
baryon rates. First, the factor (1 + 2ng(g)/n 2 ) -1 (or (1 + 2n 2 Q,^/n 2 )~ 1 ) should enter 
the formulae to eliminate the double counting of a "propagating" (anti)valon and a pair 
of the identical (anti)valons it encounters (actually, this factor may be significant for 
(anti)nucleon production only); and second, one should keep in mind that what is taken 
experimentally as A(A) is actually A(A) + £°(S) (since S° (E°) decays into A(A) + 7 in 
~ 10~ 19 s), and that the final momentum phase spaces for A and £ production differ from 
each other (compare to the above remark on n/ir 1 - '-yield ratio) as well as the combinatorial 

factors for production of the neutral and charged hyperons (see the diagram h in Fig. 2)Q 

7 The results appear to be rather insensitive to the value of (j) within this domain. 
8 As a result, the rate of A (A) production turns out to be almost twice as large as that of both 
Y,^ (E ± )-meson species. 
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Thus, we arrive at 



n 3 - 



N ~ N ~ — N- ~ N- ~ — (7) 

2 

iV A + iV s+ + AT S - ~ raQW 3 Q(g) ~ 1,6JV A (8) 

2 

iV A + iV 2+ + iV 2 _ ~ „ 1;6A r A (9) 

n z 

iV=o ~ At- ~ W&l, iv, ~ iv, + ~ ^ (io) 
2n 2 ' " " 2n 2 K ' 

N a - ^ % N« + ^ (11) 

and so on. 

The comparison of the above predictions for baryon and antibaryon rates with the 
experimental data asks for a more careful discussion (especially, for antibaryons at AGS). 
Theoretical consideration of antibaryon production is a subtle point of any model because 
it implies the details of the evolution dynamics of hadronizing nuclear matter . In contrast 
to mesons, which, once being produced in course of hadronization, preserve their identity 
(undergo predominantly elastic scattering only), the (anti)baryons start to annihilate right 
away after production; the lower is T#, the larger is the annihilation cross section. A rather 
high density of surrounding baryons (large chemical potential at AGS) and very large 
annihilation cross section, a p p ~ 100 mb, favor to very strong absorption which crucially 
diminish the antibaryon number before hadrons scatter away and stop interacting. Then, 
what remains passes through the detectors. Of course, this effect is missed in the above 
formulae. Since any theoretical estimate of this absorption could be only qualitative and is 
far from being reliable, we try to draw it directly from the experimental data. We call the 



attention to the fact [13] that measured at AGS p/ir + yield in low participant number (~ 
20) collisions is about twice as large as it is in the high participant number (~ 80-rl00) 
ones. Meanwhile, the antiproton production is a rare process and thus it is expected 
to rise up proportionally to the number of participants squared. If so, then one would 
expect that relative p/ir + yield should increase linearly. That is why one can deduce that 
annihilation eats at AGS about up to 90% of the initially produced antiprotons. At SPS, 



the relevant ratio of yields is Lki about 0,5 and thus annihilation eats (if any) not more 
than 50% of antiprotons. This is due to both substantial increase of the pionic fraction 
amongst the producing hadrons and decrease of baryonic chemical potential (both "dilute" 



the baryonic content and make lower the probability of baryon-antibaryon collisions). At 
RHIC and LHC, absorption of (anti)baryons seems getting, most probably, out of the 
game. 

Below, the results are presented in the Table (see the column th{) and Fig. 3. The 
ratios nQ(q) , WQ(s) , " Q(g) and value of (j) were found by minimization of the sum of squared 

"(5(9)' »Q(,)' n Q(s) 

deviations: 

i=k i 

X 2 * E(! - ^) 2 

i=0 "errp 

were is the number of different processes under consideration. 

Note, that till now no specific assumptions about the properties of Q7rK gas were 
made. The assumption of its quasi-ideality becomes unavoidable, when we try to extract 
the hadronization temperature Th from the ideal-gas formulae: 

^ = exp(2^ Q(g) /T), ^ = exp(2u, Q{s) /T) 

n Q(g) n Q(s) 

and 

n Q{q) _ ( m Q(g) )3/2 cx ^ ( m QW ~ m Q(g)) ] CXf3 [ ^Q(g) ~ /^QW) ] 

n s m Q(s) T ' T 

This way of doing results in nearly the same temperature, Tjj ~ (110±5) MeV for AGS, 
SPS and RHIC. Predictions for LHC given in the Table are obtained for the same value 
of T H at the point /i, = fi = 0, see Fig. 3. 
Discussion and Concluding Remarks 

It is worth mentioning that the above results for SPS are quite similar to those which 
one can carelessly obtain in the framework of the ideal gas approximation for description 
of Q7rK-phase instead using the detailed balancing equations (compare the columns thi 
and thi in the Table. This suggests that either the valon-pion-kaon gas is really quasi- 
ideal or the relative content of different components in this gas is a crude characteristic 
(this seems quite plausible) which is rather insensitive to its fine tuning. Anyway, the 
compatibility of the results seems indicative for the validity of the approach itself. 

The value of hadronization temperature, Th — 110 Mev, we have found is quite close 
to the value of thermal freeze-out temperature which has been estimated [15| for SPS and 



RHIC to be Tf ~ (100-f-110) MeV by fitting the transverse momentum spectra within 
the framework of the hydrodynamical model. If so, then fireball evolution from Th to Tf 
takes a rather short time (of some few Fermi), and one can qualitatively understand, why 
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the antiproton absorption in course and just after hadronization is crucially significant for 
AGS only (when the relative content of the nucleons is very high). In particular, that is 
why we do not need to look for a special mechanism of levelling the absorption in contrast 
to what should be necessarily invented Jl6| in the schemes with early chemical freeze-out. 

Still one point which deserves mentioning is the approximate coincidence of the hadroniza- 
tion temperatures for AGS, SPS and RHIC, although the relevant chemical potentials 
differ substantially and thus the relevant EoS's are expected to be different as well. Actu- 
ally, we can not insist on the validity of this result: \ 2 minimization is not the meaningful 
test of the theory for AGS and RHIC, since the number of parameters to be varied is 
only slightly less, than the number of the necessary data from these accelerators available 
at present. In addition, one should not take seriously too low values of \ 2 because of a 
rather insufficient accuracy of the data themselves. 

Possible variation of hadronization temperature along with valonic chemical potential 
(and thus with colliding ion energy) is closely related to the very interesting hypothesis 
that color, namely valonic, deconfmement is reasonably expected to take place already 
at very low interaction energies (ion kinetic energy being about 300 MeV only!) because 
the nucleonic density is only about twice as large as the nucleus one, and one should 
make a rather weak effort for getting valons unable to "recognize their original nucleons". 
If it is the case, then, in particular, a thermal treatment of hadron production should 
be applicable even at such a low energies. This would provide the direct indication on 
validity of the very notion of valon itself . Of course, a high-luminosity machine equipped 
with high-precision detectors is needed to try looking for the relevant manifestations. 

Thus, our conclusion is that the suggested approach provides quite successive treat- 
ment of the available data for the hadron yields as well as for the low-mass dilepton 
production || in heavy ion collisions. At the same time, it is free of the crucial inconsis- 
tency inherent in some conventional approaches, although some significant questions still 
remain which are unanswered. 

This is our pleasure to express our deep gratitude to E.L. Feinberg for the invaluable 
discussions and permanent interest to the work. 
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APPENDIX 

1. Being averaged over the particle distributions, the detailed balancing equation 
reads: 



v Q(q) (T)n„(T) ~ u^T)n Q{q) (T) (A.l) 

where Q{ are the mean values of the corresponding final state phase spaces^. Since each 
antivalon Q(q) of a certain color and flavor encounters a valon Q(q) with the probability 
n Q{g)/( n + n 7r)' the rate of Q(q)Q(q) collisions is 

***> - (1 " b) y nQW f y (a-2) 

(t) being the mean free time between the successive collisions. Quite similarly, a 7r°- 
meson encounters another 7r-meson with the probability 6, the total rate of 7T°7r collisions 
being, therefore, | fen^ || (Tr ^ 111 collisions) plus ^ J| (7r°7r° collisions); the rate of ix + n~ 
collisions is obviously ^ Of course, 7r + 7r + and ix~ix~ collisions are out of the game 
in the detailed balancing principle equation (within the above approximation), since they 
never result in a two-valonic final state. Thus, for the total rate of tttt collisions to be 
accounted one gets 

T dt 

dv„ = — bn n — — (A. 3) 

18 (t) K ' 

The averaged valonic and pionic phase spaces are 

~ Nj(2S Q + l) 2 N c pf^ and ^ ~ (2I n + l) 2 ~pj (A.4) 

respectively, where Sq is the valonic spin and 1^ is the pionic isospin, p^ is the valonic 
or pionic momentum in the CMS of two interacting valons or pions, N c and Nf are the 
color and flavor numbers (N c stands in eq. (A3) instead of N%, since only the color-singlet 
sector of the total two valon phase space is allowed for). The straightforward averaging 
over the Boltzmann distribution gives for the mean energy of a particle with mass m: 

EP(m,T) =T 2 [3^-^ + 12 + ^] (A.5) 



where Ki ;2 are the corresponding Bessel functions. The CMS value of p\ (Pq^) of each 

particle in the pionic (valonic) final state is obtained obviously by insertion into this 

9 The binary reactions are to be considered only because 2 — ► 4 reactions are substantially suppressed 
by scarcity of the typical thermal final state phase space at T ~ Th ■ 
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expression m = tuq^ (m = m n ), subtraction m 2 { m Q( q )) an d taking a half of this 
difference. Combining eq.s (Al)-(A5), we easily obtain 

[1 + 2,5(^=pL)- 1 / 2 ]- 1 (A.6) 

where the values of pf are to be taken in the accordance with eq. (A5). 

Now, consider the equilibrium kaon fraction and derive eq. (2) for the in-Q7rK-medium 
f^°' + -meson yields. For doing that, one should take into account the diagram shown in 
Fig. 1(b), 1(c). It is easy to check that, at the relevant temperatures, the averaged CMS 
momenta squared are nearly equal for all four states under consideration, and thus the 
corresponding detailed balancing equation readsQ 



(gig + 1) {2I K + 1)(2/. + l)n Q{s) n Q(q) (2I n + 1) 

n ^ + (2S Q + ly Ac " n ^ nK (A ' 7) 

(2S Q + 1) 2 N C 

n n n K 



(2I K + l)(2/ 7r + 1) 

where K stands for K + or K° and Jj is the relevant isospin. This gives eq. (2). 

Making use of these formulae, one can easily estimate that the X-meson induced 
correction to the value of b is really quite small: say, for SPS, y ~ +(4 5) 10~ 2 and it 
is undoubtedly absorbed into inaccuracy of the approach itself. 
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Figure captions 

Fiq. 1. The processes taken into account for getting the detailed balancing equations 
(see Appendix) which regulates the relative content of different components in Q-ztK gas. 
Diagrams (b) and (c) refer, actually, to the same diagram linked to the two cross channels. 

Fig. 2. The typical processes which proceed at the hadronization phase transition 
(T ~ Th). The dotted and wave lines refer to a number of the attendant pions and 
gluons, the latter being absorbed by medium. 

Fig. 3. The (yuT) phase diagram allowing for the Q7rK phase. The solid line (1) 
refers to the QGP — > QirK phase transition, the strip (2) between two dashed lines marks 
predicted in the present work temperature interval around T = 110 MeV for the expected 
value of Q7rK — > H phase transition temperature T H , and two lines (3) are quoted from 
the ref. || where they refer to the strip found by the authors for the suggested unified 
description of the hadronic chemical freeze-out after the direct QGP — > H phase transition. 
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http://arXiv.org/ps/hep-ph/0201109vl 



This figure "fig2.gif" is available in "gif" format from: 
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thi i2 


HJT H = 1.4 
// S /T ff = 0.59 
T H =H0±5 MeV 


H q /T H = 0.47 
l_i s /T H = 0.16 
T H =H0 ±5 MeV 


V q /T H = 0.026 
/i s /T H ~ 
TV=110 ±5 MeV 


fJt g /T H ~ 
^/T H ~ 
T H =H0MeV 


ThHD 




T/=110 MeV 


T /= 100 MeV 




th* 


/i g /T c/l = 1.47 
^ S /T c ^= 0.6 
T c/l =110 MeV 


fi q /T ch = 0.45 
fi s /T ch = 0.16 
T c ^=170 MeV 


fi q /T ch = 0.04 
H s /T ch < 0.01 
T c ^=190 MeV 


- 

- 



th*: model of refs. [2] and [3] ; 

thi and th?: two versions of our model, ( 7")i =0.5, 0.7 < ?9 < 1; (see in the text); 



